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A YARD PERFORMANCE MODEL BASED ON 
TASK FLOW MODELING 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional Application 
Number 60/173,436, filed December 29, 1999, which is hereby incorporated by 
reference in its entirety. 

BACKGROUND OF THE INVENTION 

This invention relates generally to railyard management, and more 
5 particularly to a yard performance model for expediting and simplifying the process of 
moving railcars through a railyard from arrival to departure. 

A terminal area is primarily used to reorganize incoming freight into 
new trains, which then move the freight to further destinations. A terminal area 
comprises one or more switchyards and interconnecting rails, and the performance of 
10 the entire terminal area depends primarily on the activities within the terminal 
switchyard(s) and the efficient dispatch of traffic within the terminal area. Thus, a 
terminal area in effect can be visualized as a small railroad network in and of itself, 
with the primary activities of managing the dispatch of traffic within the terminal area, 
and building trains within the switchyard(s). 

15 Railroad resources and activities are often analyzed hierarchically, 

from a network, regional, and local perspective. The regional perspective primarily 
emphasizes the intra-regional flow of traffic between terminal areas within a region, 
with the primary purpose of assuring that trains meet their schedules. As such, the 
regional perspective of the terminal area is not concerned with the details of train 

20 building in the terminal area, but with the overall ability of the terminals to meet the 
schedule for the flow of trains into and out of the terminal area. 

Rail corridors and terminals present an alternating sequence of services 
to railcars. The services alternate between transportation and routing. These services 
and the resources required are in limited supply. The objective of the regional level of 
25 network planning is to allocate the available resources and services within a given 
time period to optimize the flow of a collection of trains. This requires optimizing 
over time the allocation of a very large collection of factors, such as, rail segments, 
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switches, crews, locomotives, yard personnel, and yard facilities. Such an 
optimization is computationally infeasible. However, the regional view of yards and 
corridors is feasible using modeling techniques, which provide accurate, but not 
necessarily perfect, views of the capacities of yards and corridors and the relevant 
5 services involved in moving trains through the system. 

The train building process in a switchyard requires the use of tightly 
coupled and limited resources, with difficult constraints to be met as to train arrival 
and departure schedules. To model the process in detail, including the order of all 
operations, requires extensive mathematical development, and is dependent on yard 
10 topology. Therefore, it would be desirable to provide a method of allocating the 
available resources and services within a railyard in a given time period to optimize 
the flow of a collection of trains. 

The regional concept of the switchyard provides a simplified 
switchyard model posed as a sequence of car flows between reservoirs, representing 
1 5 the sequence of tasks performed on each car during the (TBP), with each reservoir 
having a limited capacity, and the flows between reservoirs being dynamically 
modulated over time by various factors present in every yard. 

BRIEF SUMMARY OF THE INVENTION 

In an exemplary embodiment, a regional terminal model focuses on the 
train building process (TBP) within a switchyard, and treats the TBP as a linear flow 

20 of cars from task to task within the yard. The model ignores inferences related to the 
specific ordering of activities in a rail yard, and assumes a first-in, first-out order of 
processing. Therefore, the primary performance metric derivable is a function of time 
that indicates how many cars can move through the yard in a specific time interval. 
Thus, a flow rate of trains in and out of the yard can be verified as possible or not 

25 possible, based on whether or not the yard can accommodate all of the incoming cars 
and process the cars at a specific time to build an outbound train scheduled to leave at 
that time. In this way, a regional overview corridor and yard use is utilized to assess 
whether or not a given yard can support a desired train schedule. More specifically, a 
regional overview is used to determine whether a given yard can support a flow of 

30 enough cars to build and depart the scheduled trains in a desired interval. 



The regional terminal model, in effect provides an envelope for yard 
capability, such that, by reordering activities within the yard, a Yard Master can 
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permute an assumed first-in, first-out order of train building, provided the total 
number of cars to be departed in the affected interval is not increased. More 
specifically, the model determines the capability of a switchyard to build trains, based 
on the schedules for arriving and departing trains and by identifying the tasks through 
5 which each car must pass in order to move from yard input to yard output, and 
modulating the flow of the cars based on yard topology and yard labor availability. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a diagram of a system used to implement the railyard 
performance model in accordance with one embodiment of the present invention; 

Figure 2 is a diagram of a railyard for illustrating the train building 
1 0 process in which the system shown in Figure 1 is utilized; 

Figure 3 illustrates the railyard shown in Figure 2 as a series of flows 
between reservoirs with valves interposed between reservoirs; 

Figure 4 shows a graphical representation depicting exemplary task 
loading in a railyard spread over 36 hours; 

15 Figure 5 shows a graphical representation depicting exemplary 

performance data for a railyard; 

Figure 6 shows flow a chart of the yard performance model; 

Figure 7 shows a flow chart of the yard performance model logic for 
train departures; 

20 Figure 8 shows a flow chart of the yard performance model logic for 

train departure yard inspection task flow; 

Figure 9 shows a flow chart of the yard performance model logic for 
classification-to-departure yard task flow; 

Figure 10 shows a flow chart of the yard performance model logic for 
25 train receive yard inspection flow; 

Figure 1 1 shows a flow chart of the yard performance model logic for 
inbound train flow; and 
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Figure 12 shows a flow chart of a sub-algorithm for sorting which of 
two choices is best for labor reassignment. 

DETAILED DESCRIPTION OF THE INVENTION 

Figure 1 is a diagram of a system 10 for implementing a yard 
performance model in accordance with one embodiment of the present invention. 
5 System 10 includes a computer 14, a display console 18 for viewing information input 
to and output from computer 14, and a user interface 22 for inputting information, 
parameters and data to computer 14. Computer 14 includes a processor 26 for 
executing all functions of computer 14, a memory storage device 30 for storing data 
and algorithms, and a database 34 for storing specific additional data. A yard master 
10 utilizes user interface 22 to input queries, parameters and data related to yard 
performance. In response to the yard master's inputs, computer 14 utilizes processor 
26, memory 30, and database 34 to solve equations and execute algorithms 
implemented in the yard performance model. 

Figure 2 is a diagram of a railyard layout for illustrating particular 
: 1 5 railyard activities involved in implementation of the yard performance model in which 
the system shown in Figure 1 is utilized. A railyard includes various sets of tracks 
dedicated to specific uses and functions. For example, an incoming train arrives in a 
receiving yard 30 and is assigned a specific receiving track. Then at some later time, 
a switch engine enters the track and moves the railcars into a classification area, or 
- . 20 bowl, 34. The tracks in classification yard 34 are likewise assigned to hold specific 
blocks of railcars being assembled for outbound trains. When a block of railcars is 
completed it is assigned to a specific track in a departure yard 38 reserved for 
assembling a specific outgoing train. When all the blocks of railcars for the departing 
train are assembled, one or more locomotives from a locomotive storage and receiving 
25 overflow yard 42 will be moved and coupled to the assembled railcars. A railyard 
also includes a service run through area 46 for servicing railcars, and a diesel shop and 
service area 50 used to service and repair locomotive. The organization of yards 
normally includes a number of throats, or bottlenecks 54, through which all cars 
involved in the train building process (TBP) must pass. Throats 54 limit the amount 
30 of parallel processing possible in a yard, and limit the rate at which the sequence of 
train building tasks may occur. Receiving yard 30, classification yard 34, storage yard 
42, departure yard 38, run through area 46 and diesel shop area 50 are 
comprehensively referred to as subyards. 
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Receiving yard 30 has a single lead-in track, which limits the rate at 
which cars can enter the TBP, and likewise for cars being moved to the receiving yard 
from overflow area 42. The number of leads between receiving yard 30 and 
classification yard 34 is likewise usually limited to one or two, as are the number of 
5 leads between classification yard 34 and departure yard 38, and the number of leads 
from departure yard 38 to a main line (not shown). Due to bottlenecks 54, there 
would normally be at most one or two switch engine crews devoted to moving 
cars/trains between the subyards. Therefore, the levels of flow for each task of the 
TBP are limited to a few options. Furthermore, use of two, rather than one, switch 
10 engines occasionally results in periods of inefficiency, due to one engine being 
blocked by the activities of the other. Therefore, the effect of multiple engines 
supporting one task will not be linear. 

However, not all tasks of the TBP directly involve engines, also 
referred to as locomotives. For example, cars arriving in receiving yard 30 must be 
1 5 inspected for defects, or "bad orders", which must be repaired before the car leaves 
yard 30. Trains being built in departure yard 38 must have all brake hoses between 
cars joined, and the brake line must be pressurized and tested before the train can 
depart. These tasks benefit linearly from additional personnel. 

Each car that arrives in the railyard to pass through the TBP requires 
20 four or five distinct tasks to be performed in serial order before it is ready as part of an 
outbound train. Each car (1) must be pulled into receiving yard 30 from the main line 
or overflow area 42, (2) must be inspected in receiving yard 30, (3) must be moved 
from receiving yard 30 to classification yard 34, (4) must be pulled, as part of a block, 
from classification yard 34 to departure yard 38, (5) must have brake hoses attached to 
25 adjoining cars and have the brakes pressure-tested as part of an outbound train. 

Figure 3 illustrates a railyard as a series of flows between reservoirs 
60, 64, 68, 72, 76 and 80 with valves 84, 88, 92, 96, 100, 104 and 108 interposed 
between reservoirs. Reservoirs 60, 64, 68, 72, 76 and 80 represent subyards and the 
valves represent modulation of the flow between subyards. Overflow or surge yard 
30 reservoir 60 is used to accommodate incoming trains only if there is insufficient space 
in the receiving yard reservoir 64. Receiving yard reservoirs 64 and 68, and departure 
yard reservoirs 76 and 80 are each depicted, in Figure 3, as two reservoirs, 
dichotomized according to whether railcars have or have not been inspected. The 
receiving yard is for bad orders and railcars that have or have not been inspected in 
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the departure yard to verify that the brake hoses have been coupled and pressure- 
tested. Thus, receiving yard reservoir 64 represents railcars in receiving yard 30 
(shown in Figure 2) that have not been inspected, and receiving yard reservoir 68 
represents railcars in receiving yard 30 that have been inspected. Likewise, departure 
5 yard reservoir 76 represents railcars in departure yard 38 (shown in Figure 2) that 
have not been inspected, and departure yard reservoir 80 represents railcars in 
departure yard 38 that have been inspected. 

Each of valves 84, 88, 92, 96, 100, 104 and 108 denote a task flow 
modulation associated with a modulating agent. Inbound flow valve 84 is effectively 

10 opened for inflow by the T-Plan (PLEASE EXPLAIN T-PLAN), and surge yard 
reservoir 60 must be prepared to accept the increase in level. Surge-to-receiving flow 
valve 88 is modulated by the need to move railcars from surge yard reservoir 60, and 
the availability of an engine crew to effect that action. Receiving inspection (RI) flow 
valve 92 is modulated by the availability of a carman to perform inspections and a 

1 5 hostler for removing power from incoming trains. Classification or bowl flow valve 
96 is modulated by engine crews and a pin puller, who are actively moving railcars 
from receiving yard reservoir 68 to bowl reservoir 72. Departure flow valve 100 is 
modulated by longfielder(s) and engine crew(s). Departure inspection (DI) flow valve 
104 is modulated by brakemen and hostler(s) who couple and inspect brakes and 

20 attach power to the trains. Finally, outbound flow valve 108 is modulated by the T- 
Plan departure schedule for the yard. 

In addition to the direct modulation of the flows by labor and the T- 
Plan, the regional yard model includes flow rate modulations related to (1) internal 
yard congestion that affects engine movements when a subyard is nearly full, (2) 
25 dynamic reassignment of labor by a yard master, in order to alleviate the more severe 
backlogs, and (3) shift-dependent modulation of labor rate during a shift, which is 
commonly seen in yards. Each of these modulating influences will be discussed 
below in detail. 

Other primary factors, which affect the output of a railyard over a 
30 specific time window are (1) initial conditions in the yard, such as the backlogs of cars 
in each yard, the proportion of receiving and departure cars already inspected, (2) the 
total amount and distribution of labor available for a shift, (3) maximum possible flow 
rates for all of the tasks, for example, a receiving yard with two lead-in tracks may be 
able to accept more cars per hour than one with one lead-in track, and (4) the nominal 
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rate at which each task will be carried out, given nominal crewing. Parameters for all 
these inputs is gathered or inferred from yard archival data. 

Although activities in a switching yard are continuous, the yard 
performance model advances the state of the yard in discrete time increments, for 
5 example, an interval At will pass during which flow rates will remain constant and 
subyard car levels will vary accordingly. A yardmaster may then examine the subyard 
levels and modulate the flow rates according to a labor policy. The interval At will be 
kept short, for example, 15 minutes, so that yardmaster decisions occur before any 
large imbalances build in the railyard. 

10 Although the mechanics of the yard performance model are simple, 

very substantial oscillatory behavior can develop related to the level of cars in each 
subyard, the flow rates between subyards, the processes by which the flow is 
modulated, and the amount of delay involved in recognizing an imbalance and taking 
corrective action. Short-term behavior of a railyard will have an unpredictable and 

1 5 high variance, while long term overall flow through the yard is relatively predictable 
with a much lower variance. 

Utilizing the yard performance model, an analyst enters relevant 
parameters, simulates yard task flow through a specified time interval, and determines 
if the train schedules for the yard during the associated time window can be met. 

20 Yard Parameters 

Table 1, below, shows a list of the data parameters, and exemplary 
data, associated with the yard performance model. 



TABLE 1 
Table of Yard and Train Parameters 



Scenario 

Yard parameter file: 
Input schedule file: 
Output schedule file: 
Simulation length: 
Initial clock time: 
Analysis window: 
Random seed: 
Labor algorithm: 



North Platte W 

InSched 

OutSched 

36.0hr. 

8:00 

0:00-8:00 
12345 

dynamic, backlog first 
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Scheduling: 
Congestive effect: 
TDM: 



fixed 
active 
active 



Yard Topology and Initial State 

Surge Receiving Bowl Departure 

Capacity (cars): 250 1000 600 1200 

Initial Levels: 500 300 600 

% Inspected: 80 80 

Congestion effect - Starts at 90%, slows by 20% 

Train Arrival Parameters 

Max. pull-in speed (mph) : 15.0 

Mean car length (ft.): 60.0 

Mean interarri val time (hr . ) : 1.3 

Interarrival std. dev. (nr.): 0.25 

Mean train length (cars): 1 06.0 

Train length std. dev.: 22.0 



Labor Parameters 

Arrivals to yard: 
Surge to receive yard: 
Receive yard inspection: 
Classification: 
Bowl to departure yard: 



Max. cars/hr. 
660 
300 
300 
240 
600 



Departure inspection & power: 1000 
Labor 



Actual Cars/hr per schedule 



60 
66 
79 
72 



Mix 


Surge 


Receive 


Classify 


Trim 


Departure 


EC 


0.0 


0.0 


1.0 


1.0 


0.0 


HS 


0.0 


0.0 


0.5 


0.0 


0.5 


CM 


0.0 


1.0 


0.0 


0.0 


0.0 


PP 


0.0 


0.0 


1.0 


0.0 


0.0 


LF 


0.0 


0.0 


0.0 


1.0 


0.0 


BM 


0.0 


0.0 


0.0 


0.0 


1.0 



Shift start: 7:00, TDM offset: 6.0 hrs., TDM decrease: 20% 



A yard parameter file, an input schedule file, and an output schedule 
file are file names associated with the railyard and train schedules under analysis, and 
are prestored in system 10, or defined on-screen during a session. Once defined, the 
45 yard parameter, input schedule and output schedule can be saved within system 10 for 
later reuse. The simulated length and initial clock time entries are based on the local 
yard clock time and the intended length of time for which yard task flow is to be 
simulated. The analysis window parameter specifies a period of time, referenced to 
clock time, during the simulation when incoming cars will be counted. A time 
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interval during which those cars exit the yard is shown in output graphs of the 
program. The random seed parameter is used to initialize a random number generator 
and uses various modes of operation. The labor algorithm parameter reflects a yard 
master policy for dealing with congestion in the TBP tasks. The algorithm selected 
5 may be "static", indicating no labor movements during the simulation, "dynamic, 
headend first", indicating that backlogs at the front of the task flow are given priority 
over backlogs toward the end, or "dynamic, backlog first", indicating that higher 
backlogs will take priority over lower backlogs. The scheduling parameter is either 
fixed, denoting that a specific, defined schedule is in use, or random, indicating that 

10 the analysis program can generate random schedules for arriving and/or departing 
trains. The congestive effect parameter relates to when a subyard is nearly full so the 
engines involved in the tasks for that yard may be slowed by the need to choose a less 
preferred route between points. The congestive effect parameter is either active or 
inactive during the analysis. The TDM parameter refers to a time-dependent 

1 5 modulation of task rates, which occurs in many yards. The TDM is the effect of a 
general slowdown in labor rate toward the end of a shift. This is specifiable in the 
program in terms of an offset from the beginning of the shift when the effect is noted, 
and a percentage by which the task rates slow. 

The yard topology parameters are based on the capacities of the 
20 subyards, in terms of railcars. The initial levels of cars in each yard are determined 
before any simulation begins based on those cars in the receive and departure yard. A 
percent of railcars will already have been inspected (receive yard) or ready for 
departure (departure yard) at the time that simulation begins. 

The train arrival parameters provide a maximum pull-in speed 
25 parameter, which reflects the upper speed limit at which an arriving train may enter 
the receiving yard. The maximum pull-in parameter affects the maximum number of 
cars that can enter the receiving yard in any specified interval. The mean car length 
parameter, is entered in feet. Along with the maximum pull-in speed, the mean car 
length determines how many cars can enter the yard in a specific interval. The mean 
30 interarrival time paramter and interarrival standard deviation parameter are used if 
random scheduling has been chosen. If random scheduling is chosen, the train 
arrivals and departures will be normally distributed according to the chosen values for 
mean and standard deviation. The mean train length parameter and train length 
standard deviations parameters are used to assign normally distributed train lengths, 
35 given in cars, to trains which are generated by the random scheduling process. 
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The labor parameters reflect the rates at which tasks can be done in the 
yard, and the mix of crew members nominally assigned to each of the five flow tasks. 
As shown in Table 1, there are maximum possible rates and actual rates. The 
maximum rates are functions of the physical topology of the yard, and reflect 
5 limitations not of labor, but of available tracks. The actual rates are the task flow rates 
for each task as a function of crews actually assigned to the task at the initiation of the 
simulation. Typically, the actual flow rates for the tasks will never approach the 
maximum limits. 

The labor mix provides a matrix denoting how labor is initially 
10 assigned to the tasks. The labor mix may vary during simulation if a dynamic labor 
assignment algorithms has been chosen. The abbreviations of labor categories used in 
the matrix are as follows. 

EC — yard engine crew participating in hump and trim tasks. 
HS - hostler, relating to removing power from the incoming trains and 
15 delivering power to outbound trains. 

CM -carman inspecting for bad orders in the receiving yard. 
PP - pin puller, relating to decoupling cars for the classify process. 
LF - longfielder, relating to assisting the classification process by 
correcting anomalies in the classification yard. 
20 BM- brakeman working in the departure yard to couple brake hoses 

and pressure test a completed train. 

The labor parameters are user-specified during program execution. 
However, there are a number of restrictions on labor assignments. For example, labor 
25 is not completely fungible, such that, the hostlers, engine crews, and the pin pullers 
cannot change job categories. However, the carmen, brakemen, and longfielders can 
be freely reassigned within those three categories. 

Finally, there are three other parameters that describe the TDM 
process. Time-dependent modulation of task rates is represented as a decrease in 
30 labor rate at a certain time after the beginning of the shift. Thus, the time of shift 
initiation, relative to the clock time for initiation of the simulation, the time at which 
task flow rates decrease, as an offset from the beginning of the shift, and the percent 
by which the task flow decreases from the offset to the end of the shift, are necessary 
parameters. 

35 In addition to the parameters described above, schedule parameters 

defining a train schedule need to be defined. Schedule parameter are typically saved 
to and retrieved from a disk for use in the yard performance mode. In one 
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embodiment, the pertinent schedule parameters are the number of trains the schedule 
will contain, the total time interval for the schedule, whether or not inbound train 
arrival times will be perturbed normally about their nominal values, a number 
assigned as a train ID, the number of cars in the train, the expected arrival time of the 
5 train, offset from the beginning of the simulation interval, and a standard deviation 
about inbound train arrival time, which will be used with a normal distribution to vary 
train arrival times about their specified arrival times when the variable input schedule 
mode is chosen. 

The standard deviation about inbound train arrival is not the same as 
10 the standard deviation of interarrival times or the standard deviation of train length. 
In the case of standard deviation of interarrival times and train length, entirely random 
schedules are generated corresponding to the means and standard deviations of train 
interarrival times and train lengths. In one embodiment, a defined schedule is used, 
but the arrival times of the trains may vary somewhat around the specified arrival 
1 5 times. Thus, a train arrival time could be perturbed to be greater than the arrival time 
of the next train in a sequence. When trains stay in sequence, the yard performance 
model slides later trains out when a given train arrival time will impact later trains. 
However, it only delays later trains as required, not all later trains. 

Also, if a schedule time is shorter than the specified simulation 
20 interval, the program appends randomly generated arrivals and/or departures to the 
schedule in order to continue a flow into and out of the yard. 

Furthermore, other parameters, all of which are not specified by the 
user, are utilized in the yard performance model. The backlog level at which 
congestive effects begin for any subyard and the maximum decay of task rate as 

25 congestion moves to 100% are parameters specified by the user. Since yard 
congestion typically does not affect the task rate of a carman, pin puller, longfielder, 
or brakeman, backlog level and maximum decay only apply to the proportion of each 
task rate which is contributed by engine crews. If an engine crew operates between 
two congested yards, for example, pulling blocks from classification to departure, then 

30 the congestion factor from both yards affect the associated task rate. The minimum 
backlog at which a yardmaster will consider moving labor, the borrowing limit below 
which a task backlog must be before the yard master can move labor from the task, 
and the task flow rates, are not user defined parameters. The borrowing limit 
parameter is typically set to 75% and only applies to the dynamic, headend-first labor 
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assignment algorithm. For dynamic, backlog-first, labor can be borrowed from any 
task with a lower backlog than the borrowing task. Task flow rates are functions of 
the labor assigned to each task. Equations internal to the program provide this 
functional relationship and the forms and coefficients of the equations are not 
5 accessible to the user. 

Figure 4 shows a graph representation 200 depicting the task loading in 
a yard spread over 36 hours, based on the exemplary initial conditions shown in Table 
1, above. A receiving loading bar graph 204, a classification loading bar graph 208, 
and a departure loading bar graph 212 depict subyard loading over time. The bars 

10 indicate backlogs below 60%, between 60% and 80%, between 80% and 100% , and 
at 100% saturation. The sawtooth effect in both the receiving yard bars and the 
departure yard bars reflects the discrete moments at which trains arrive and depart, 
respectively. Small squares 216, superimposed on receiving loading bar graph 204 
and departure loading bar graph 212, reflect the percentage, as a fraction of each bar, 

15 of the cars that have been inspected for bad orders in receiving yard 30 (shown in 
Figure 2), or are ready for departure from departure yard 38 (shown in Figure 2). A 
labor management graph 220 depicts the labor management process for the duration 
of the simulation. At the left of graph 220 is an axis subdivided into equal intervals to 
represent each task. The abbreviations shown in graph 220 are given as, 

20 S-R - moving cars from the surge yard to the receiving yard, 

R&I - inspecting and removing power in the receiving yard, 
R-C - moving cars from receiving into the classification yard, 
C-D - moving blocks of cars from classification to departure, 
P&I - coupling brakes, adding power, and pressure-testing in departure. 

25 

The discontinuous bars of graph 220, shown for each task, represent 
labor assigned to each task over time. The varying thicknesses of the bars in graph 
220 indicate more or less labor of the corresponding category assigned to a task. 
Graphs 204, 208, 212 and 220 provide a detailed look at the yard for the period of 
30 simulation, which aids in the understanding of what yard activities may be having a 
deleterious effect on throughput. 

Figure 5 shows a graphical representation 300 depicting exemplary 
performance data for a railyard in terms of its ability to service incoming and outgoing 
trains as per desired regional train schedules. A yard output graph 304 shows a yard 
35 output versus time plot, representing the number of cars departed from the yard over 
the time of the simulation. Graph 304 is effectively a step function because the cars 



-12- 



GEH-0 1-067 



are departed in discrete segments corresponding to the trains in an output schedule. A 
graph 308 indicates train arrivals and departures across the time interval of the 
simulation. The lengths of the vertical bars of graph 308 are proportional to the length 
of the train, with full scale corresponding to the maximum number of cars that could 
5 enter/leave a railyard in the simulation update interval The full scale is based on the 
mean pull-in speed and car length. For graph 308, full scale is 165 cars. In graph 
308, the bars above the horizontal axis represent arriving trains, and the bars below 
the axis represent departures. A train departure cannot occur on schedule if departure 
yard 38 (shown in Figure 2) does not contain enough inspected cars at the scheduled 
1 0 departure time of the train. In that case, the train will be delayed. 

Parametric Description of the Flow Model 

Figure 6 shows flow chart 350 of the yard performance model. 
(PLEASE VERIFY CHANGES TO FIGURE 6 ARE CORRECT IN 
ACCORDANCE WITH THIS PARAGRAPH AND THE NEXT) The model is 

15 best described by partitioning the performance model into submodels. Once initial 
conditions, such as train schedules 354, initial backlogs 358, yard topology 362 and 
labor assignment 366 are input, the model calculates 370 the initial task flow rates 
based on an initial state as input by the user. A user, such as a yard master, utilizes 
user interface 22 and display console 18 (shown in Figure 1) to access to all 

20 parameters of the model, except the non-user specified parameters discussed above, 
and may modify the default parameters either by editing during program execution, or 
by recalling previously saved files for train schedules and yard parameters. T 0 , of the 
initial state is the clock time at the yard in which the simulation begins. Next, the 
model updates 374 task backlog of each of the five tasks discussed in relation to 

25 Figure 2 above, and computes or modifies 376 task flow rates. For example, the 
model advances cars to the next task, based on the flow rates in effect. This process 
begins with train departures, and works backward to the beginning of the yard. The 
task flow rates are updated in accordance with the varying yard conditions. Task 
backlog updates and task flow rate updates are done on a time increment of fifteen 

30 minutes, so that each task moves a corresponding number of cars to the next task. 
After each task backlog is updated, flow rates are updated, according to one or more 
flow modulating effects, such as, modifying 378 engine crew task rates, modifying 
382 all task rates, and activating 386 a new labor mix. At the end of each task flow 
rate update, the time is checked 390, and, if it equals the end time of the simulation, 
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the update loop ceases 394 and outputs the graphics shown in Figures 4 and 5 to 
display console 18 (shown in Figure 1). 

If the simulation time has not reached the endpoint, the model 
determines 398 if the time is an even hour, and if either of the dynamic labor modes 
5 have been chosen. If so, the yard master labor decision process is executed 402, 
whereby labor may be moved between tasks according to the labor assignment 
restrictions discussed above. However, the decision to move labor incurs 406 a 
fifteen minute penalty, during which time the reassigned personnel are in transit to the 
new assignment. If the time-dependent modulation mode has been selected 410, then 

1 0 the simulation time is compared to shift time to determine if the task flow rates should 
be modified, for example, decreased or increased. If any of the subyards have become 
congested 414 since the last update the appropriate adjustments are made to engine 
crew task rates affecting the engine crews working in the affected subyards. After the 
flow rate adjustments have been made, the simulation clock time is updated 420 by 

15 advancing the time by 15 minutes, and the main loop beginning with updating 374 
task backlogs is repeated. 

Referring to Figures 2 and 6, in one embodiment, updating task 
backlogs 374 involves updating the backlog of five tasks. For example, each 
incoming car which will be placed in an outbound train (1) must be pulled into 
20 receiving yard 30 from the main line or the surge yard (not shown), (2) must be 
inspected in receiving yard 30, (3) must be moved from receiving yard 30 to 
classification yard 34, (4) must be pulled, as part of a block, from classification yard 
34 to departure yard 38, and (5) must have brake hoses attached to adjoining cars and 
have the brakes pressure-tested as part of an outbound train. 

25 The yard performance model handles the task flow for these five tasks 

in reverse order, for example, the first task is to depart a train if the schedule so 
warrants, updating the backlog in the departure yard, and then the program works 
backward through tasks and backlogs to the first task. 

Referring to Figure 3, a railyard is modeled as comprising six separate 
30 subyards, with both receiving yard 30 and departure yard 34 (shown in Figure 2) 
effectively comprising two subyards, one of cars inspected and a second of cars not 
inspected. The cars inspected in the receiving yard are referred to as cars in the RI 
subyard, and similarly, cars already inspected in the departure yard are referred to as 
cars in the DI yard. Although there are six logical subyards, there are seven flows, 
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corresponding to the seven valves 84, 88, 92, 96, 100, 104 and 108, to modulate. 
These flow valves are denoted as inbound flow valve 84 , surge-to-receiving flow 
valve 88, RI flow valve 92, classification flow valve 96, departure flow valve 100, DI 
flow valve 104, and outbound flow valve 108. Inbound flow valve and outbound 
5 flow valve are not implicitly controlled by the yard, but represent the inbound and 
outbound trains scheduled by explicit predefined schedules. However, valves 88, 92, 
96, 100 and 104 represent tasks that are modulated directly under the control of the 
yard. 

In reference to Figures 7 through 17, the following variables are 
1 0 defined in accordance with the following list. 

(A) Let Cj, C2, C3, C4, C5, C(5 be the maximum capacities, in cars, of 
each of the six subyards, for example, the surge, receiving, RI, 
classification, departure, and DI subyards, respectively. 

(B) At any time t, let Lj(t), L2(t), 13(f), Ltft), L$(t), L$(t) be the 
1 5 occupancy, in cars, of each of the six subyards, 

(C) Let i?/, R2, R3, R4, R5, R6, Rj represent the maximum flow rates 
possible for each of the six subyards. 

(D) At any time t, let F](t), F 2 (t), F 3 (t), F 4 (t), F 5 (t), F 6 (t), F 7 (t) be 
the actual flow rates, in cars per minute, for each of the six subyards, 

20 and let Dt represent the time increment (15 minutes) for each yard 

update. 

(E) Let A j, A2, Ana, be the trains scheduled to arrive during the 
simulation period, with T^(Aj) being the arrival time for train Af. 

(F) Let D], D 2 , ... , Dnd, be the trains scheduled to depart during the 
25 simulation period, with Tj)(Dj) being the scheduled departure time for 

train Dp 

(G) Let L(H) denote the length of any train H, in cars. 

(H) Let Eft) represent the modulation factor applied to task flow rates 
at time t corresponding to TDM. 

30 (I) Let N2, N3, N4, N$, N$ be the nominal flow rates for the internal 

yard tasks (cars/minute), given nominal crewing, for the six subyards. 
(J) Let Uj(t), U 2 (t), U3(t), U4(t), Us(t) be the congestion factors for 
each of the engine crew or hostler time t, where the congestion factor is 
a value in the range [0,1] which indicates slowdown of engine crew or 

35 hostler activity due to the subyards being nearly or completely full. 

More specifically, 

U](t) = slowdown of an engine crew for the surge-to- 
receive task, 

40 U 2 (t) = slowdown of the hostler for the receive-to- 

classify task, 

U3(t) = slowdown of the engine crew for the receive-to- 
classify task, 

U4(t) = slowdown of the engine crew for the classify-to- 
45 departure task, and 

U5(t) = slowdown of the hostler for the departure yard 
inspection task. 
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Figure 7 is a flow chart 450 of the yard performance model logic for 
train departures. When utilizing the yard performance model, a train departure 
depends on the train departure being scheduled at the corresponding moment of 
5 simulation time, and there being sufficient inspected cars in the departure yard to 
make up a train of the required length. T], T2,..., Tj, Tj+ ], represent the sequence of 
discrete times at which all flows are updated. The yard performance model logic for 
train departures first checks 454 the departing train schedule to see if a train is due to 
depart at time Tj. If not, no further train departure logic is executed, and the model 

10 advances 456 to departure inspection task logic. Otherwise, train Dj is scheduled 458 
for departure. Next it is determined 462 if the number of cars L(Dj) is less than or 
equal to the number of cars L$(Tj) ready for departure in the DI subyard. If no, then 
there are sufficient cars for the departing train, and it is departed, resulting in 
decrementing 466 the DI subyard level by the value L(Di). If the number of cars 

1 5 L(Dj) is not less than or equal to the number of cars L$(Tj) ready for departure in the 
DI subyard, train Dj is delayed 470 by a time set 474 to at least time Tj+ j, and any 
following trains which might be departing in consecutive time slots from Tj+j 
outward will be pushed back one slot. The model then advances 456 to departure 
inspection task flow logic. 

20 Figure 8 is a flow chart 500 of the yard performance model logic for 

train departure yard inspection flow. The departure yard train inspection moves trains 
from departure yard reservoir 76 to the DI yard reservoir 80 (shown in Figure 3) 
according to the current flow rate for the inspection task, modified 504 by the TDM 
factor. The only constraint is that the total backlog of inspected cars cannot exceed 

25 508 the total number of cars in departure yard 38 (shown in figure 2). If so, the model 
sets 512 total backlog of inspected cars equal to the total number of cars in departure 
yard 38, and the model advances 516 to classification-to-departure task flow logic. 

Figure 9 is a flow chart 520 of the yard performance model logic for 
classification-to-departure yard flow. The basic classification-to-departure yard car 

30 flow attempts to move the number N of cars corresponding 524 to the flow rate and 
simulation update time interval At. However, N is subject to two constraints. N can 
not exceed 528 the number of cars available in classification yard reservoir 72 (shown 
in Figure 3). If N exceeds the number of available cars, the model sets 532 N equal to 
the number of cars available in classification yard reservoir 72. Additionally, 

35 departure yard reservoir 76 (shown in Figure 3) must have adequate room to 
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accommodate 536 N more cars. If the departure yard reservoir 76 does not have 
adequate room, the model sets 540 N equal to the available room. N is then modified 
544 downward, if necessary, the car levels for classification reservoir 72 and departure 
yard reservoir 76 are then updated accordingly, and the model advances 546 to 
5 receive-to-classification task flow logic. 

The basic receive-to-classification yard flow is similar to the 
classification-to-departure yard task flow shown in Figure 9. Receive-to-classification 
flow attempts to move the N cars corresponding to the flow rate and simulation update 
time interval At, from RI reservoir 68 (shown in Figure 3) to classification yard 

10 reservoir 72. However, N is subject to two constraints. N can not exceed the number 
of cars available in RI reservoir 68, and classification yard reservoir 72 must have 
adequate room to accommodate N more cars. N is then modified downward, if 
necessary, and the car levels for RI reservoir 68 and classification yard reservoir 72 
are updated accordingly. The logic diagram for this process is identical to that of 

15 Figure 10, except for decrementing all subscripts by 1. 

Figure 10 is a flow chart 550 of the yard performance model logic for 
train receive yard inspection flow. The receive yard train inspection moves trains 
from receive yard reservoir 64 to the RI yard reservoir 68 (shown in Figure 3) 
according to the current flow rate for the inspection task, modified 554 by the TDM 
20 factor. The only constraint is that the total backlog of inspected cars cannot exceed 
558 the total number of cars in receiving yard 38 (shown in figure 2). If so, the model 
sets 562 total backlog of inspected cars equal to the total number of cars in receiving 
yard 38, and the model advances 566 to surge-to-receive task flow logic. 

The logic for the surge-to-receive yard car movement task is identical 
25 to that for the classification-to-departure yard tasks discussed in reference to Figure 
10. The surge-to-receive model logic is demonstrated by uniformly decrementing all 
subscripts in Figure 10 by 3. 

Figure 1 1 is a flow chart 600 of the yard performance model logic for 
inbound train flow. The yard performance model logic for inbound trains first checks 
30 454 the inbound train schedule to see if a train is due at time Tj. If not, no further 
inbound train logic is executed, and the model advances 608 to the task flow update 
logic (shown in block 374 of Figure 6). Otherwise the receive or surge yard car count 
is updated 612 by the number of cars in the arriving train. Typically, an inbound train 
goes to receive yard reservoir 64 (shown in Figure 3) if there is room, otherwise the 
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inbound train goes to surge yard reservoir 60 (shown in Figure 3) if there is room, and 
failing either of those alternatives, the train is not accounted for, for example, it is 
regarded as left out on the main line. Next it is determined 616 if the number of 
inbound cars is greater than the capacity of receiving yard reservoir 64. If not, then 
5 there is sufficient room, and the train is moved to receiving yard reservoir 64, the 
level of receiving yard reservoir 64 is incremented 620, and the model advances 608 
to task flow update logic. If the number of incoming cars is greater than the capacity 
of receiving yard reservoir 64 the model determines if the number of incoming cars is 
greater than the capacity of surge yard reservoir 60. If not, then there is sufficient 
10 room, and the train is moved to surge yard reservoir 60, the level of surge yard 
reservoir 60 is incremented 628 and model advances 608 to task flow update logic. 

In Figure 6, blocks 402 and 386 indicate that a yard master may 
reallocate labor on every hour, measured from the beginning of simulation time, if 
either of the dynamic labor assignment modes, for example, dynamic, headend first, 
1 5 or dynamic, backlog first, has been chosen. Otherwise, the process is effectively null, 
and labor remains constant for the entire simulation. Two major factors effect labor 
allocation, (1) task flow rates as a function of labor assignment and (2) yard master 
labor allocation algorithms. 

Task flow rates as a function of labor assignment involve the yard tasks 
20 of surge-to-receive movement, receive yard car inspection, receive-to-classification 
car movement, classification-to-departure car movement, departure yard train 
building. The arrivals and departures of trains from the yard are not considered a yard 
task. Each task has a default nominal flow rate based on a nominal crewing 
assignment, and these nominal flow rates can be altered by the user to reflect 
25 conditions at any yard of interest. Associated with each nominal flow rate is a 
nominal crew assignment, which is fixed by the model. These nominal crew 
assignments are as follows, 

(1) Surge-to-Receive Car Flow ~ 1 engine crew; 

(2) Receive Yard Car Inspection — 1 carman; 

30 (3) Receive-to-Classification Car Flow — 1 engine crew, 1 

pinpuller, X A hostler; 

(4) Classify-to -Departure Car Flow — 1 engine crew, 1 

longfielder; 

(5) Departure Train Building ~ V 2 hostler, 1 brakeman. 

35 

The yard performance model expresses each task in terms of task flow 
rate equations. Each flow rate equation is expressed in terms of the nominal flow 
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rates, congestive factors, and the actual labor assignment to the task at any moment; 
the TDM modulation also affects task rates, but is handled extrinsically to the task 
flow equations. 

The surge-to-receive flow rate involves the labor category of engine 
5 crews (EC). The equation providing the flow rate for cars from the surge yard to the 
receive yard is given by 

F2(N E M=N 2 U 2 (t)f}l'-' , (1) 

10 where Nj?c(t) — the number of engine crews assigned to surge-to-receive car flow at 
time t. 

The summation of powers of two in this equation provides a 
diminishing return for additional engine crews. The first crew assigned counts as 1 
crew, the next as Y 2 of a crew, and so forth, because multiple engine crews will be 
1 5 somewhat self-congestive. This self-congestive effect will be used in some of the 
following equations. 

When considering the receiving yard car inspection flow rate, the 
carman operating in the receiving yard are not regarded as affected by a yard 
20 congestion factor, or by the presence of more carmen operating in the receiving yard. 
Therefore, the equation for the inspection task flow rate is linear, and given by 

F 3 (N CM (t)) = N 3 N CM (t), (2) 

25 where Ncm(0 = me number of carmen assigned to receive yard car inspection at time 
t. 

The receiving-to-classification yard car flow depends on the presence 
of a pinpuller, an engine crew, and the hostler who is removing the power from 
incoming trains. Without an engine crew or pinpuller, there can be no flow at all into 
30 the classification yard. However, the absence of the hostler decreases the flow, but 
does not force it to zero, because the road crews arriving with the inbound train may 
be induced to move power from the incoming trains to the service area. Reflecting 
these facts, the equation for car flow is given by 
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if N EC if) = 0 or N pp {t) = 0, then flow rate is zero, otherwise 



f n f ( Neck) 2N HS (t) > 

F< (N EC {t\ N PP (/), N HS (/)) = I °J 0.4 + 1 0.4C/ 3 (f ) £ 2'" 1 + 0.2J7 2 (f ) £ 2'" 1 



(3) 



where N]?c(t) = number of engine crews assigned to classification at time t, Npp(t) = 
number of pinpullers active at time t, and NhsO) = number of hostlers operating at the 
5 receive end of the yard at time t. 

The classification-to-departure yard flow requires an engine crew and, 
optionally, a longfielder. The absence of the engine crew results in a zero flow rate 
even if a longfielder is present. The equation for this flow is given by 

ifN]7c(t) = 0 then flow rate is zero, otherwise 

io ^(^cW^W)^ 0 ^^^,)^-.^^^ , (4) 

where Ngcft) = number of engine crews assigned to the trim process at time t, and 
NlfCO = number of longfielders assigned to the trim process at time t. In this case 
the nominal labor rate is realized with just an engine crew assigned, and the 
1 5 longfielder increases the nominal car flow rate by 20%. 

The departure yard train building flow depends on brakemen, to couple 
brake hoses and pressure test the brakes, and a hostler to mate power to the departing 
trains. In the departure yard train building flow, the brakeman is indispensable. 
Therefore, the flow rate is zero if there is no brakeman assigned. However, the hostler 
20 is not absolutely necessary, since the road crews assigned to the departing trains may 
also bring power to the train. The equation for departure yard train building flow is 
given by 

{0 If NBM(t) = 0 , then the flow rate is zero, otherwise 



25 where NfiMft) = number of brakemen assigned at time t, and Nffs(t) = number of 
hostlers assigned at time /. 
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The yard master labor allocation algorithms of the yard performance 
model define the manner by which labor is reassigned during a simulation. The 
model permits three choices of labor allocation, static, headend first, and backlog first. 
In the static model labor is not reassigned at all during the simulation. 

If labor is reassigned, it is reassigned to alleviate backlogs in the tasks, 
because such backlogs can lead to congestive effects in the yards, or to a task being 
unable to move cars forward through the process due to the next task having no 
storage room remaining for incoming cars. Therefore, headend first or backlog first 
labor allocation is used. Headend first method gives priority to backlogs near the 
front of the yard before backlogs toward the end of the yard. The reasoning behind 
this method is that if the front of the yard is saturated, then inbound trains remain out 
on the mainlines or sidings, which is the most undesirable situation. 

The backlog first method addresses the maximum backlog first, where 
backlog is defined as the percent by which a task is backlogged relative to the total 
capacity of cars that the corresponding task can accommodate. 

Since inspection tasks differ from the actual car flow tasks in terms of 
what backlog means, the following equations express the backlogs at time / for each 
task. The backlog for the surge-to-receive task is given by 

B 2 (t) = L/fl/C,. (6) 
The backlog for the receive yard car inspection task is given by 

B(t) = l° if ^(0 = 0 

3W [1 - Z ? (/•)/ Z, 2 (/) otherwise ' ^ ; 

The backlog for the receive-to-classification car flow backlog is given by 

B 4 (t) = L 2 (t)/C 2 . (8) 
The backlog for the classification to departure yard car flow is given by 

B 5 (t) = L 3 (t)/C 3 . (9) 
The backlog for the departure yard train building task is given by 



-21- 



GEH-0 1-067 



*H-4vW <10) 

The above definitions of backlogs are used for either of the dynamic 
labor assignment algorithms, for example, the headend first method or the backlog 
first method. The dynamic labor assignments uses the following definitions, 

(1) Bi = the backlog level at or above which an attempt will be made 
to obtain labor for the corresponding task (B£ = 0.8), and 

(2) Bg = the backlog limit at or below which it is acceptable to borrow 
labor from the corresponding task (Bg = 0.75). 

Labor categories within the yard are not completely fungible. The 
matrix in Table 2, below, indicates what labor categories are interchangeable, and 
from which tasks labor of a given category may be drawn when a reassignment 
occurs. The headings across the top of the matrix represent the five tasks of moving 
cars the surge yard to the receiving yard, moving cars from the receiving yard to the 
RI yard, moving cars from the RI yart to the classification yard, moving cars from the 
classification yard to the departure yard, and moving cars from the departure yard to 
the DI yard. The labor available in any cell may be reassigned to the other tasks with 
an 'X' in the same row. 



Table 2 



lasks 





S-R 


RI : R-C 


C-D 


Dl 


EC 


X 


X 


X 




HS 






X 




X 


CM 




X 




X 


X 


PP 






X 






LF 




X 




X 


X 


BM 




X 




X 


X 



As shown in Table 2, when labor is to be reassigned, it is, in all cases 
except one (hostlers), obtainable from exactly two other tasks. Therefore, the model 
algorithm that reassigns labor uses a subalgorithm, which sorts out which of two 
choices is best for the reassignment. 

Figure 12 is a flow chart 650 of the subalgorithm for sorting out which 
of two choices is best for labor reassignment. The subalgorithm accepts 654 as input 
-22- 



GEH-0 1-067 



two labor categories, two tasks, and the backlog of the task for which labor is being 
potentially recruited. Next the subalgorithm determines if either or both of the tasks 
from which labor may be recruited have any labor of the required category 658, or 
have backlogs less than a specified backlog 662. Both of these criteria must be met 
5 for a task to be able to supply labor 664. If both criteria are not met the subalgorithm 
sets 668 the task choice backlog to a value of 10. Next the task choices are compared 
672 and the one with the lower backlog is chosen 676. If neither of the task choices 
have labor that may be recruited or have backlogs less than a specified backlog, the 
subalgorithm returns 674 a negative task number, signifying that labor of the 
10 categories desired cannot be borrowed from the candidate tasks. The subalgorithm is 
denoted by 

L=(T 1t T 2 ,L 1r L 2 ,B) (11) 
where Tj and Tj? are the tasks from which labor might be borrowed, L ] and L2 are 
15 the labor types to be borrowed, respectively from Tj and T2, and B is a backlog level 
below which a donor task can permit borrowing. 

Also, the indexing of tasks and labor categories in the following discussion will 
correspond to the entries in the matrix of Figure 13, with the upper left (blank) cell 
corresponding to (0,0). 

20 

The headend first algorithm (HFA) checks the backlogs of the tasks, in 
order from the front of the yard, for example, surge-to receive car movement, to the 
end of the yard, for example, brake coupling and inspection in the departure yard. 
The HFA then attempts to move labor from other tasks to aid any task with a backlog 
25 greater than B£ = 80%. However, the attempted borrowing process differs somewhat 
from one task to the other. Borrowing of labor will only be permitted from tasks with 
backlogs below Bq = 75%. 

As shown in Table 2, enhancing the flow rate from surge-to-receiving 
yard requires additional engine crews. Thus, labor in the form of engine crews (labor 

30 type 1) is borrowed either from the receive-to -classification task or the classification 
to departure task. The subalgoritm function is called to determine which (if either) of 
the tasks might supply the engine crew, returning that choice as the value of the 
variable "min". If min < 0, no labor can be obtained. Otherwise, the affected labor 
matrix entries LM(ij) are modified to reflect movement of an engine crew from the 

35 donor task to the surge-to-receive yard task. 
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The receive inspection task can borrow either a longfielder or a 
brakeman, converting those labor categories to the carman category. The 
subalgorithm determines which, if either, of the longfielder or brakeman will be 
borrowed, and, if either is available, the labor matrix is updated accordingly. 

5 The receive-to-classification yard flow differs from the previous cases 

in that it normally relies on two labor categories, engine crews and hostlers. The 
engine crews can be obtained from either of two other tasks. Thus, the subalgoritm 
looks for an engine crew in surge-to -receive task or the classify-to-departure task, and, 
if a crew is available, updates the labor matrix accordingly. 

10 The search for additional support from the hostler is different, since the 

hostler(s) are only involved in one other task, and the subalgorithm function does not 
apply. However, typically a yard has one hostler acting to support both tasks 
involving hostlers. Thus hostlers can be moved in fractions of l A, for example, the 
single hostler serving two tasks may be placed exclusively on one task, changing his 

1 5 effect in increments of 0.5 . 

The classification-to-departure yard car flow depends on longfielders 
and engine crews. This task flow has the property that it will be zero (no car flow) if 
there is no engine crew, even if a longfielder is available. Therefore, the logic for this 
process is to first obtain, if possible, an engine crew, and then only look for 
20 longfielder support if the task has at least one engine crew assigned. A longfielder is 
sought only if at least one engine crew will be assigned to this task, and that 
longfielder can be obtained by converting a carman or brakeman to longfielder duty. 
Thus, the subalgorithm function is applied to determine which, if either, category of 
labor will be converted. 

25 Additionally, in the classification-to-departure task, a longfielder 

without an engine crew has no value and the task flow is zero. Therefore a yard 
master should not remove the last engine crew from the classification-to-departure 
task without reassigning the longfielder as well. 

The departure inspection task depends on both hostlers and brakemen. 
30 A brakeman can be obtained by converting either a longfielder or carman to brakeman 
service. Therefore, the subalgorithm function is used to choose from where the 
additional brakeman is acquired. 
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As a final action of the HFA, the subalgorithm checks to see if the 
longfielder is without an engine crew, and if so, the longfielder is converted to a 
carman or brakeman. The specific alternative chosen is to favor the recipient task 
with the highest current backlog. 

5 The only distinction between the HFA and the Backlog First Algorithm 

(BFA) is that the HFA proceeds through the tasks in car flow order to determine 
whether to move labor, while the BFA first reorders the tasks so that labor is 
reassigned to tasks with the order of reassignment decision corresponding to the task 
backlogs in descending order. Thus the subalgoritym for the BFA first creates an 
1 0 ordered list of the tasks, corresponding to descending backlogs, and then applies the 
same logic as applied in the HFA, except that the order of labor reassignments follows 
the ordered task list. 

However, there is one other distinction. The HFA uses fixed backlog 
thresholds to determine which tasks may be donors, and which recipients, 
15 corresponding to the backlog levels Bg = 75% and Bi = 80%, respectively. The BFA 
uses the same threshold for recipient tasks, but allows a donor task to be any task with 
a backlog lower than the current backlog of the recipient task. 

The labor reassignment in either HFA or BFA cannot be cyclic, such 
that, one task may not borrow and then later be made to return a crew member within 
20 the same reassignment cycle. The ordering of donor versus recipient backlog 
thresholds prevents this from occurring. 

There is a delayed effect of labor reassignment in that a labor 
reassignment decision pays a penalty of transport time for those crew members 
moving to new tasks. Thus, the new values of the labor mix do not take effect 
25 immediately, but are delayed by one simulation cycle, of At =15 minutes. Therefore, 
in the next iteration of task flow after the labor reassignments, the labor contingent 
will be lower than the entire crew available in the yard, which will be reflected in the 
flow rates for that time interval. Thereafter, until the next labor reassignment, the 
flow rates will reflect all yard labor, as reassigned. 

30 The congestive effect only applies to engine crews and hostlers, and 

only affects task rates in proportion to the contribution of the engine crew or hostler to 
the overall task rate. The congestion factors applicable to the five yard tasks 
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involving engine crews and hostlers, are included in the task flow rate equations (1) 
through (5). 

The congestive effect is determined by two parameters, given as 

Cg = a level of task backlog at which the congestive effect first begins, 
and 

Cj = the maximum decay in engine crew/hostler efficiency caused by 
congestion, when the corresponding subyard backlog reaches 
subyard capacity. 



10 Each engine crew operates between two subyards, and if both of those 

subyards are congested, the congestive effects of both subyards are applied as 
multiples to the nominal engine crew contribution to the task. Thus, the values 
Uj(t),... Us(t) are themselves products of two subyard congestion factors. Although 
the entire yard is divided into six subyards in the flow model, the RI and DI subyards 

15 do not create congestion effects, since the backlogs in these yards are not directly 
related to the physical capacity for car storage in any physical subyard. For each of 
the physical subyards, let 

F](Lj(t)/C]) = congestion factor for the surge yard at time t, 
F2(L2(t)/C2) = congestion factor for the receiving yard at time t, 
20 F4(L,4(t)/C4) = congestion factor for the classification yard at time t, 

and 

F5(Ls(t)/C5) = congestion factor for the surge yard at time t. 
(PLEASE CLARIFY IF THE LAST EQUATION SHOULD 

HAVE SUBSCRIPTS OF '6' AND SAY "FACTOR FOR THE DEPARTURE 

25 YARD") 

These factors are parameterized in terms of the their current backlogs, 
and each factor is defined by the equation 



F,Mt)/C,) = 



1 if Z,(/)/C, <C B 

'AAwdf ' (12) 



30 where (for /= 1,2,4, 5). 

Therefore, the congestive effect has no effect until the threshold level 
Q? is reached, and then descends as a cosine curve from 1 to the maximum effect of 
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1 - Cjr over the remaining part of the backlog range. Additionally, the factors 
Fj(Lf(t)/Ci) must be related to the values Utft). Specifically, 



The TDM Effect 

The time dependent modulation (TDM) effect, represents a variation in 



apparent labor rates during the course of a shift. It is an optional effect in the yard 
performance model. The TDM effect represents the fact that normally, the labor rates 
are higher in the early part of a shift, and decrease toward the end of the shift. It is 
parameterized in the program in terms of three parameters, 



Tg - time of shift initiation for any shift, 
Tj) = time after beginning of shift that labor rates decrease, 
p = the proportional labor rate decrease from the labor rate observed 
in the early part of the shift. 

Using a nominal task flow rate of r for any task, for example, an 



average over a full shift, the task flow rate for any task must be elevated above the 
nominal rate assigned by a certain amount for the period of time before the rate 
decreases, so that the overall mean task rate does conform to the nominal rate listed on 
the parameters page. Specifically, if R hi represents the higher than nominal rate 
observed in the early part of a shift, and R lo represents the lower than nominal rate 
observed in a latter part of the shift, then for a shift of length T L , 




(13) 




(14) 



which can be expressed as 




from which 
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T D (l-p)+T L p jand 

(15). 

R,o=PK 

are obtained. 

In an exemplary embodiment, using the TDM parameters specified in 
Table 1, a nominal task rate r will be elevated to 1.0526r for 6 hours, and then will 
decrease to 0. 8421 r for 2 hours. The mean task rate is then [6(1. 05r) + 2(.84r)]/8 = 
r, and .84r/1.05r = 0.8. During each task flow update in the simulation, the yard 
performance model determines whether the current simulation time is in the high task 
rate, or the low task rate part of a shift. If the simulation time is in the early part of 
the shift, all task rates are elevated from the nominal rates corresponding to r, as 
defined for equations 14-15, to Rfa as given by equation 15 and for the latter part of 
the shift, all task rates are decreased from their nominal rates r to i?/ 0 , as given by 
equation 15. 

While the invention has been described in terms of various specific 
embodiments, those skilled in the art will recognize that the invention can be practiced 
with modification within the spirit and scope of the claims. 
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